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ABSTRACT 

Although collisionless shocks are ubiquitous in astrophysics, certain key aspects of them are not well understood. 
In particular, the process known as collisionless electron heating, whereby electrons are rapidly energized at the 
shock front, is one of the main open issues in shock physics. Here, we present the first clear evidence for efficient 
collisionless electron heating at the reverse shock of Tycho's supernova remnant (SNR), revealed by Fe K diagnostics 
using high-quality X-ray data obtained by the Suzaku satellite. We detect K/i (3p — > lx) fluorescence emission 
from low-ionization Fe ejecta excited by energetic thermal electrons at the reverse shock front, which peaks at a 
smaller radius than Fe Ka (2 p -+ lx) emission dominated by a relatively highly ionized component. Comparisons 
with our hydrodynamical simulations imply instantaneous electron heating to a temperature 1000 times higher 
than expected from Coulomb collisions alone. The unique environment of the reverse shock, which is propagating 
with a high Mach number into rarefied ejecta with a low magnetic field strength, puts strong constraints on the 
physical mechanism responsible for this heating and favors a cross-shock potential created by charge deflection 
at the shock front. Our sensitive observation also reveals that the reverse shock radius of this SNR is about 10% 
smaller than the previous measurement using the Fe Ka morphology from the Chandra observations. Since strong 
Fe K/i fluorescence is expected only from low-ionization plasma where Fe ions still have many 3 p electrons, this 
feature is key to diagnosing the plasma state and distribution of the immediate postshock ejecta in a young SNR. 

Key words: atomic data - hydrodynamics - ISM: individual objects (SN 1572, Tycho's SNR) - ISM: supernova 
remnants - shock waves - X-rays: ISM 
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1. INTRODUCTION 

A supersonic flow colliding with another medium will form 
a shock wave, which can be described as a discontinuity in the 
physical conditions of the flow. Shock waves have been exten- 
sively observed on Earth and also in a number of astrophysical 

settings with a wide range of spatial scales: from the solar wind 
(e.g., Tidman & Krall 1971; Schwartz et al. 1988) to afterglows 
of gamma-ray bursts (GRBs; e.g., Meszaros & Rees 1997; Sari 
et al. 1998) and merging galaxy clusters (e.g., Markevitch et al. 
2005; Markevitch & Vikhlinin 2007). Unlike terrestrial shocks, 
the low densities found in some astrophysical environments 
imply that the shock transition occurs on scales much shorter 
than the typical particle mean free path for Coulomb collisions. 
Therefore, collisionless processes involving collective interac- 
tions between particles and electromagnetic fields must be re- 

sponsible for the shock formation (McKee 1974). Despite the 
ubiquity and importance of these collisionless shocks in as- 
trophysics, the detailed physical processes that determine their 
fundamental properties are still poorly understood. 

One particularly mysterious process, which might be closely 
related to the formation of the collisionless shock front itself, 
is the rapid heating of electrons at the shock transition. A 

strong shock with velocity v s will result in a downstream 
temperature 7} = 3 in; i> 2 /16A:b, where m, is the mass of 
particle i and k K is the Boltzmann constant. Since the timescale 
for collisional equilibration between different species is much 


longer than the time a particle spends in the shock transition zone 
(Spitzer 1962), the equation above can be applied independently 
to each species i, unless collisionless processes contribute 
to the temperature equilibration. It follows that the electron 
temperature (T e ) should be much lower than the temperature 
of heavier ions (7] on ) immediately behind the shock and they 
will slowly equilibrate to a common temperature via Coulomb 
collisions further downstream. However, a number of theoretical 
investigations have suggested that rapid collisionless electron 
heating can occur at shock fronts (e.g., McKee 1974; Cargill & 
Papadopoulos 1988; Laming 2000; Ghavamian et al. 2007). 

Supernova remnants (SNRs) offer an ideal site to study 
this heating observationally, because they form long-lived, 
fast shocks in both the interstellar medium (ISM) and the 
supernova ejecta. To date, most work has concentrated on 
Balmer-dominated shocks associated with SNR blast waves 
expanding into the ISM (e.g., Raymond et al. 1983; Laming 
et al. 1996; Ghavamian et al. 2001, 2007; Rakowski et al. 2003; 
Helder et al. 201 1). Studies of electron heating in reverse shocks 
(RSs), which are propagating into the supernova ejecta are, on 
the other hand, very limited (Hamilton et al. 1997; France et al. 
2011). In this paper, we focus on the RS of Tycho's SNR, the 
remnant of the Type la supernova observed in A.D. 1572. 

There are a number of important differences in environment 
between blast waves and RSs. Unlike the ISM, supernova ejecta 
generally have only a small fraction of neutral particles, because 
the strong ultraviolet/soft X-ray flux from the shocked material 
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Figure 1. Suzaku XIS images of Tycho's, SNR in the (a) 6.43-6.53 keV (Fe Ka) 
band and (b) 7. 7-9.0 keV band (continuum emission). North is up and east is to 
the left. The northwest (NW) region confined with the green lines is where we 
extract the spectrum shown in Figures 2 and 4. 

(A color version of this figure is available in the online journal.) 

can easily photoionize the heavy elements of the unshocked 
ejecta (e.g., Hamilton et al. 1997). Furthermore, in a Type la 
SNR the ejecta consist only of heavy elements, with no con- 
tribution of hydrogen and helium. The magnetic field strength 
expected in the interior of the white dwarf before the explo- 
sion ( B < 10 13 G; Suh & Mathews 2000) will be dramatically 
diluted by expansion. From flux conservation, a conservative up- 
per limit on the field strength present near the RS location in Ty- 
cho's SNR is estimated to be ~ 1 0“ 7 G. This is sufficiently high 
that the ion gyroradius, r g ~ 10 13 (A/56) (Z? / 10 7 [G]) _1 cm 
(where A is the mass number), is much smaller than the SNR 
radius of ~10 19 cm, ensuring the formation of a collisionless 
RS. Yet, the inferred strength is at least an order of magnitude 
below the typical value in the ISM (a few times 10~ 6 G). The 
geometry of the interior magnetic field should be highly ordered 
as the field is stretched radially by the expanding ejecta. These 
distinct properties, strongly contrasting with those in the ISM, 
allow us to probe the physics of collisionless electron heating 
under conditions far different from those of any earlier work. 

Previous hydrodynamical calculations applied to Tycho's 
SNR required the presence of a modest amount of collisionless 
electron heating at the RS to explain the observed X-ray flux 
from the shocked ejecta (Badenes et al. 2005, 2006). This is 
in contrast to earlier results on another Type la SNR, SN 1006, 
where little evidence for collisionless electron heating was found 
(Hamilton et al. 1997). Since the RS in Tycho's SNR has begun 
to propagate into the Fe-dominated ejecta (e.g., Hwang et al. 
1998; Warren et al. 2005; Badenes et al. 2006), the emission lines 
from Fe directly probe the conditions (i.e., electron temperature) 
in the postshock region. In general, the initial shock-heated 
material in SNRs is at a very low charge state and only gradually 
becomes ionized by collisions with hot free electrons. The ejecta 
in young SNRs are, therefore, in a non-equilibrium ionization 
(NEI) condition (Masai 1984). Collisional interactions between 
hot electrons and low-ionized Fe in an NEI plasma produce inner 
K-shell ionization of the Fe ions, followed by Ka(2p —> Is) 
or K/3 (3 p Is) fluorescence transitions (e.g., Palmeri et al. 

2003; Mendoza et al. 2004). These emission lines are excellent 
diagnostics of the temperature and charge population in the 
shocked material. 

Here, we present new, strong evidence for collisionless 
electron heating at the RS of Tycho's SNR, revealed by 
Fe K-shell X-ray lines from sensitive observations with the 
X-ray Imaging Spectrometer (XIS) on board the Suzaku satellite. 
Our observations allow the first-ever detailed study of weak Fe 
K/3 emission alongside the stronger Fe Ka line. In Section 2, we 



Figure 2. XIS spectrum of Tycho's SNR in the 5. 0-9.0 keV band, where the 
background data are subtracted. A phenomenological model using a power law 
(dashed line for the continuum) and five Gaussians (solid gray lines for the 
emission lines) yields the best-fit parameters given in Table 1 . 

analyze observational data based on the state-of-the-art atomic 
physics models we compute. In Section 3, we constrain the ef- 
ficiency of collisionless electron heating by comparing our hy- 
drodynamical calculations and finally we discuss the plausible 
mechanism of this efficient heating process. 

2. OBSERVATIONAL RESULTS 

The observations with the Suzaku XIS were performed during 
2008 August with a total effective exposure of 415 ks. The 
primary data reduction was performed in accordance with the 
standard procedure recommended by the instrument team. We 
use only the data of front-illuminated CCDs (XIS0 and 3) and 
merge them to improve the photon statistics. Figure 1 shows XIS 
images in the energy bands of the Fe Ka line ((a) 6.43-6.53 keV) 
and continuum emission ((b) 7. 7-9.0 keV). This continuum has 
been shown to be dominated by synchrotron radiation from 
relativistic electrons accelerated at the blast wave (Hwang et al. 
2002; Cassam-Chena'i et al. 2007). We extract an X-ray spectrum 
from the northwest (NW) rim indicated in Figure 1 . Background 
data are taken from the nearby sky and subtracted from the 
source spectrum. Figure 2 is the resultant spectrum in the 
5. 0-9.0 keV band, where emission from Cr Ka, Mn Ka, Fe 
Ka and K/3, and Ni Ka is clearly resolved. We note that this is 
the first detection of Ni from this remnant, although we do not 
focus on this element in this paper. 

2.1. Fe K Emission Diagnostics 

We measure the centroid energies of the Fe Ka ( E K „ ) and 
K/3 ( E K p) line blends and their flux ratio CR) using Gaus- 
sian functions and obtain E Ku = 6435 ± 1 (±6) eV, E K /i = 
7104 ± 10 (±7) eV, and K = 5.5+° 0 6 5 %. The parenthetical values 
are the instrumental systematic uncertainty in the energy scale, 
0.1% of the mean energy, due to the incomplete gain calibra- 
tion of the XIS (Ozawa et al. 2009). The continuum is simply 
modeled by a power law, giving a photon index of T ~ 2.9. 
The best-fit parameters for all the emission lines are given in 
Table 1. During the analysis, an absorption column density of 
JVh = 7 x 10 21 cm -2 (Cassam-Chena'i et al. 2007) with a so- 
lar elemental composition (Wilms et al. 2000) was assumed, 
although the analyzed energy band was not affected by the fore- 
ground extinction. To evaluate the systematic uncertainty due to 
background subtraction, we perform fits using a different back- 
ground data set that includes only the instrumental component 
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Table 2 

Theoretical Values of the Fe Kcr and K/l Centroid Energies and their Intensity 
Ratios for the Different Charge Numbers z 


log(n e t [crrr 3 s]) ® ? 

1 J/ | i i i i | i 


10 11 

t — I— i — I — I I I I 



Figure 3. Expected centroid energies of the (a) Fe Kg' and (b) Fe K/? emission, 
and (c) the K/3 /Ka flux ratio as a function of the charge number z of Fe ions, 
with the corresponding ionization ages ( n e t ) indicated at the top. The best-fit 
values for Tycho's SNR are shown with the solid green lines. The dashed green 
lines in panel (c) indicate the lcr lower and upper limits of the observed value. 
The black squares and red circles are the values calculated using the atomic data 
of Palmeri et al. (2003) and Mendoza et al. (2004), respectively. We also use the 
FAC to calculate the rate coefficients of collisional ionization and excitation. 
The blue triangles are obtained by full calculations using the FAC (K. A. Eriksen 
et al., in preparation). 

(A color version of this figure is available in the online journal.) 


Table 1 

Best-fit Spectral Parameters for the NW Rim 


Emission 

Centroid 

FWHM 

Flux 


(eV) 

(eV) 

(10 -6 photons cm -2 s -1 ) 

Cr Key 

5482/ 1 1 ° 2 

141+35 

5 05 +0-72 
- 0.69 

MnKd 

6012/1 

141 

1 70+0.47 

1 • 7 —0.46 

FeKcr 

6435 ± 1 

138 ±2 

107 ± 1 

Fe K/l 

7104 ± 10 

160 ±42 

5 62 +0 ’ 61 
J - oz — 0.56 

Ni Key 

7478 ± 32 

138 

1 O' - ) +0.42 
0.43 


Notes. The uncertainties are the statistical component in the 1 o confidence 
range. The Gaussian widths (FWHM) of the Mn Kg' and Ni Kg' lines are linked 
to those of Cr Kg? and Fe Ka, respectively. 


(non-X-ray background). No significant change is found in the 
measurement of the line centroids and intensities. 

The observed Fe-line parameters, E^ a , E K p, and 1Z, are 
compared in Figure 3 with the theoretically expected values 
for different charge numbers z (where z = 1 indicates singly 
ionized Fe), which are also tabulated in Table 2. For z = 0-7 
and z = 8-16, we use level energies, Einstein A-values, and 
fluorescence yields provided in the archival database of Palmeri 
et al. (2003) and Mendoza et al. (2004), respectively. Since the 
transition probabilities of forbidden processes (e.g., 2s -> Is) 
are negligible in multiple-electron ions (Palmeri et al. 2003; 
Mendoza et al. 2004), we take into account only 2p Is 
and 3 p —> Is transitions as radiation channels for the Ka 
and K/l emission. We calculate rate coefficients for collisional 
ionization and excitation for each charge number using the 
“Flexible Atomic Code” (FAC; Gu 2008). For z ^ 16, we 
perform full computational calculations with the FAC to obtain 


Z 

Ena 

(eV) 

Enp 

(eV) 

n 

z 

Ekc 

(eV) 

Ekp 

(eV) 

n 

0 

6402 

7059 

0.120 

12 

6414 

7141 

0.069 

1 

6402 

7060 

0.121 

13 

6419 

7153 

0.052 

2 

6402 

7060 

0.122 

14 

6425 

7159 

0.022 

3 

6401 

7059 

0.127 

15 

6428 

7176 

0.010 

4 

6400 

7063 

0.132 

16 

6427 

7192 

0.012 

5 

6399 

7070 

0.136 

17 

6455 

7270 

0.009 

6 

6399 

7075 

0.141 

18 

6484 

7351 

0.013 

7 

6399 

7081 

0.149 

19 

6517 

7434 

0.025 

8 

6398 

7090 

0.168 

20 

6544 

7517 

0.029 

9 

6401 

7102 

0.146 

21 

6575 

7610 

0.036 

10 

6405 

7115 

0.122 

22 

6589 

7705 

0.044 

11 

6410 

7128 

0.096 

23 

6641 

7777 

0.075 


the theoretical values (K. A. Eriksen et al., in preparation). 
During the calculations, we assume an electron temperature of 
5 keV. The temperature dependence is found to be significant 
only for z ^ 20, where the inner K-shell excitation rate 
becomes dominant over the K-shell ionization rate. Since the 
population of such highly charged Fe is not substantial in Tycho’s 
SNR, as discussed below, the diagnostics we perform here are 
essentially independent of the electron temperature of the Fe 
ejecta. 

As found in Figure 3(a), the observed E Ka value corresponds 
to the charge states Fe l5+ ~Fe l7+ and an ionization age ( n e t ) 
of ~1 x 10 10 cm -3 s, where n e and t are the electron density 
and the time elapsed since shock heating, respectively. This 
result is consistent with several previous measurements (e.g., 
Hwang et al. 1998; Hayato et al. 2010). It has, therefore, 
frequently been assumed that all the shocked ejecta responsible 
for the Fe K-shell emission have this ionization age. We 
find in Figure 3(b), however, that the observed C K /i value 
corresponds to significantly lower charge states Fe s+ ~Fe 10+ , 
with an ionization age of ~1 x 10 9 cm" 3 s, about 10 times 
lower than that indicated by the Fe Ka centroid. The K/f-to-K a 
flux ratio (7 Z) is also sensitive to the charge number, especially 
in the range z — 8-14 (Figure 3(c)). In this regime, the flux ratio 
experiences a drastic decrease because the Fe ions lose their 3 p- 
shell electrons (which are responsible for the K/l fluorescence) 
as z increases. The observed value is closest to the expected 
ratio for z = 13, intermediate between the results from the Ka 
and K/l centroids. 

The inconsistency among the three diagnostics indicates the 
presence of a range of plasma conditions, with the Ka emission 
being dominated by more highly ionized Fe and the K/l emis- 
sion by less ionized Fe. We re-fit the NW spectrum applying 
a ‘two-component’ model for the Fe emission. The red Gaus- 
sians in Figure 4 correspond to the higher ionization component, 
where E K /j and 7 Z are fixed to 7200 eV and 1% (the values theo- 
retically expected for z ~ 16). Only the Ka centroid is allowed 
to vary, yielding = 6447 7 3 eV, which is in between the 
values for Fe 16+ and Fe 17+ . The contribution of the lower ioniza- 
tion component is indicated by the green Gaussians in Figure 4, 
where we fix E Ka and TZ to 6400 eV and 15% (as expected for 
z ~ 8). The best-fit E K p value of 7090± 1 1 eV is consistent with 
that for Fe 8+ . To summarize, we are able to explain all the Fe 
K-shell emission self-consistently with this simple two- 
component model using different ionization states and the 
expected K/l /Ka emissivity ratios. 
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Energy [keV] 


Figure 4. Same spectral data as in Figure 2, but the “two-component” model 
is applied to the Fe K emission (see the text for details). Red and green 
represent the higher-ionization (around Fe ls+ ) and lower-ionization (around 
Fe 8+ ) components, respectively. The K/3 / Ka flux ratio is expected to be lower 
in the former (~0. 01) than the latter (~0.15), because Fe 16+ ions have no bound 
electron in the 3 p shell, which is responsible for the K/i emission. 

(A color version of this figure is available in the online journal.) 



Figure 5. XIS image of Tycho’s SNR in the 7. 0-7. 2 keV (Fe K/> ) band, where the 
Fe Ka? image (same as Figure 1(a)) is overplotted in contours. The synchrotron 
continuum flux estimated using Figure 1(b) is subtracted from the raw Fe K/3 
image. The morphology is subject to some uncertainties since a spatially uniform 
photon index is assumed in this subtraction procedure. It is nevertheless clear 
that the Fe K/> emission peaks at a smaller radius than the Ka emission in the 
bright NW rim. The green sector indicates where the spatially resolved spectral 
analysis (Figure 6) is performed. 

(A color version of this figure is available in the online journal.) 

2.2. Spatial Analysis 

The interpretation in the previous subsection predicts that the 
Fe K/3 emission peaks interior to the Ka emission, because the 
innermost ejecta were heated by the RS more recently and so 
should have a lower ionization age than the outer ejecta. We 
perform a spatial analysis to verify that this is indeed the case. 
Since the Fe K/3 emission is not as strong as the synchrotron 
continuum flux in the same energy band, subtraction of the 
continuum component is necessary. We estimate the continuum 
level in the 7. 0-7. 2 keV (Fe K/3) band by scaling the synchrotron 
X-ray image (Figure 1(b)) using a photon index of 2.9 (the 
best-fit value for the NW rim spectrum). The color image in 
Figure 5 is created by subtracting these scaled data from the 
raw 7.0-7.2-keV image, where contours of the Ka emission 
(corresponding to the image in Figure 1(a)) are overlaid. As we 
expected, the K/3 emission has a smaller peak radius than the 
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Figure 6. Radial profiles for the surface brightness of the Fe Ka and K/3 emission 
in the NW quadrant, derived from the spatially resolved spectral analysis. The 
smaller peak radius of the Fe K/3 is confirmed. The RS positions determined by 
Warren et al. (2005) and this work are indicated with the dashed lines. 

(A color version of this figure is available in the online journal.) 

other. It should be noted, however, that this imaging analysis has 
some uncertainties, because we assume that the photon index 
of the continuum X-rays is spatially uniform over the entire 
region. This is actually not the case, but the synchrotron emission 
in Tycho's SNR is known to have spatially inhomogeneous 
hardness (Cassam-Chenai et al. 2007; Eriksen et al. 2011). 

This motivates our spatially resolved spectral analysis, which 
gives a more quantitative measurement of the Fe K/3 distribution. 
The azimuthal sector shown in Figure 5 is divided into 10 
radial zones of variable width (to account for the variation 
in brightness) for spectral extraction, by assuming the same 
center position for the SNR as that determined by the previous 
Chandra observations (Warren et al. 2005). The spectrum from 
each region is fit with the same model applied to the NW rim 
(five Gaussians plus a power law), allowing the photon index 
as well as the other parameters to vary freely. The resulting 
surface brightness profiles of the Fe Ka and K/3 emission are 
shown in Figure 6, confirming the smaller peak radius of the K/3 
emission, while the Ka peak radius is consistent with that from 
the Chandra data. 

Our analysis demonstrates that the radial profile of the Fe K/3 
emission is more sensitive to the immediate postshock ejecta 
than the Ka emission, revealing a ~10% smaller RS radius 
than was previously determined using the Fe Ka morphology 
from the Chandra observations (Warren et al. 2005). Assuming 
a simple shell geometry, which predicts that the maximum 
surface brightness will coincide with the inner edge of the 
shell (Warren et al. 2005), we estimate the RS radius to be 
7.1 x 10 18 (D/3.0 [kpc]) cm, or 2.3 (D/3.0 [kpc]) pc, where/) 
is the distance to Tycho's SNR. This is about 63% of the SNR 
blast wave radius, 3.6 (/)/ 3.0 [kpc]) pc. 

3. INTERPRETATION 

3.1. Comparison with Hydrodynamical Calculations 

Our new X-ray measurements have revealed that Fe in the 
innermost ejecta is in an extremely low ionization state (z < 8). 
Yet, the strong Fe K/3 emission requires the electron temperature 
near the RS front be high enough to ionize the inner K-shell 
electrons of these low-z ions. Keeping this in mind, we constrain 
the efficiency of collisionless electron heating by comparing 
our results with one-dimensional hydrodynamical simulations 
that incorporate an NEI calculation (Badenes et al. 2006). For 
the initial conditions, we assume the ejecta structure expected 
for a typical delayed-detonation Type la supernova with an Fe 
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Figure 7. Electron temperature as a function of the mean charge of Fe ions from 
our hydrodynamical simulations. The corresponding radius is also given above. 
The black curve is the ftmin model where no collisionless electron heating 
is assumed. The temperature ratio between the electrons and ions at the RS 
front is, therefore, set by their mass ratio. The models represented by the red, 
blue, and green curves assume that collisionless electron heating occurs at the 
RS, parameterized by (/ 3 = T e /T l0n ) with values set to 0.003, 0.01, and 0.03, 
respectively. 

(A color version of this figure is available in the online journal.) 

yield of ~0.8 M G and an explosion energy of ~ 1.2 x 10 51 erg. 
This model reproduces the fundamental properties of Tycho’s 
supernova: the historical light curve (Ruiz-Lapuente 2004) 
and modern light echo spectrum (Rest et al. 2008; Krause 
et al. 2008). The SNR evolution is followed to the age of 
Tycho’s SNR assuming a uniform ambient density (pam) of 
2 x 10“ 24 g cm -3 . The result is also in good agreement with the 
basic dynamics (e.g., angular size, shock speed) of the SNR for 
a reasonable value of the distance (Badenes et al. 2006). The 
effect of collisionless electron heating is introduced to our NEI 
calculations by parameterizing the electron-to-ion temperature 
ratio just behind the RS, /i = T e /T lon . 

Figure 7 shows the relationship among charge state, radius, 
and electron temperature for several values of ft. The black 
curve is derived under the assumption that the initial tempera- 
ture of each species follows the equation 7) = 3 m,- u 2 /16£b, 
which leads to T f . / 7p e = m e /mp e ~ 10“ 5 at the RS front (here- 
after, the y6 m ; n model). The other models, shown as the red, blue, 
and green curves, have larger values of /3, which result in higher 
electron temperatures in the postshock region (hereafter, the col- 
lisionless heating models). Subsequent temperature changes are 
due to collisional processes, heating via ion-electron Coulomb 
collisions (dominant in the /3 m i n model) and cooling via the 
ongoing collisional ionization process (prominent in the colli- 
sionless heating models). The /-i n , m model predicts an electron 
temperature of <1 keV in the region dominated by z ^ 10 
Fe ions. The free electrons in this region are, therefore, not 
energetic enough to produce significant K-shell ionization and 
subsequent fluorescence. This is shown more quantitatively in 
Figure 8(a), where we plot the expected Fe Ka and K/l lumi- 
nosity for each charge state. The p mm model clearly fails to 
reproduce the strong Kfi emission from low-ionized Fe, in di- 
rect conflict with our observations. In contrast, the collisionless 
heating models can reproduce the emission from the broad Fe 
ion population, as illustrated in Figure 8(b) (corresponding to 
P = 0 . 01 ). 



0 5 10 15 20 

Charge Number 

Figure 8. (a) Predicted luminosities of the Fe Kck (red) and K/3 (green) emission 
lines for the /3 m i n model. Owing to the low electron temperature in the innermost 
region, little or no emission from low-ionized (z < 10) Fe is expected, (b) Same 
as panel (a), but for the model with (3 = 0.01. Emission from Fe with various 
charge states, including z < 10 where higher 7 Z values are achieved, is expected. 
(A color version of this figure is available in the online journal.) 



(3min 0.003 0.01 0.03 

P (Te/Tion at RS) 


Figure 9. Comparison between the predicted values for (a) E k<*, (b) Ej^p, and 
(c) 7 Z and the observed values, for which statistical and systematic uncertainties 
are indicated with dark and light red regions, respectively. The thin solid (with 
squares), thick solid (with circles), and dashed (with triangles) lines correspond 
to ambient densities (pam) of 1 x 10 -24 , 2 x 10 -24 , and 3 x 10 -24 g cm -3 , 
respectively. We find that the model with pam = 2 x 10~ 24 g cm -3 and 
f3 = 0.01 provides the best match to the observations. 

(A color version of this figure is available in the online journal.) 

In Figure 9, we compare the model-predicted values of E k«, 
Ej^p, and 1Z with the observed values for Tycho ’s SNR, where 
the impact of different ambient density values, pam = 1 x 10 -24 


5 


The Astrophysical Journal, 780:136 (6pp), 2014 January 10 


Yamaguchi et al. 


and 3 x 1 O ' 24 g cm -3 , are also explored. We confirm that the 
fi m in models cannot reproduce the low and high TZ values 
we observe for any value of the ambient density. Although a 
lower pam value allows for a larger population of low-ionized 
Fe, it also leads to a Fe K« centroid energy far lower than 
the observed value. Only the collisionless heating model with 
f = 0.01 and Pam = 2 x 10 -24 g cm -3 explains all the observed 
values within the given uncertainties, indicating that when the 
electrons pass through the RS front, they gain an internal energy 
about 10 3 times higher than expected without collisionless 
heating. 

It should be noted that the recent infrared observation with 
Spitzer determined the mean ISM density around Tycho's SNR 
to be (2-4) xl0~ 25 g cm -3 (Williams et al. 2013). This is 
consistent with the previous estimates from the blast wave 
expansion rate (Katsuda et al. 2010) and the upper limit of 
thermal X-ray emission from the shocked ISM (Cassam-Chena'i 
et al. 2007), but is significantly lower than the value from our 
diagnostics. This discrepancy implies the presence of a density 
gradient around the progenitor; the SNR blast wave had initially 
interacted with relatively high-density matter, which enhanced 
the ionization age of the outer ejecta but is now expanding 
into the low-density ISM. A similar interpretation was given by 
Dwarkadas & Chevalier (1998) and Chiotellis et al. (2013). 
Although the effect of such a non-uniform ISM should be 
involved in our future calculations, we believe that this does 
not affect our conclusion significantly. The Fe K emission from 
the broad range of ionization ages can be achieved only when the 
collisionless heating efficiency is high enough, as demonstrated 
in Figures 8 and 9. 

3.2. Origin of Collisionless Electron Heating 

We have presented the first clear evidence for efficient 
collisionless heating of electrons at the RS of Tycho's SNR, 
which is propagating into metal-rich (no-hydrogen) ejecta at 
low magnetic field strength, with a velocity of ~4000 km s -1 
(Badenes et al. 2006; Hayato et al. 2010). As the unshocked 
ejecta in young Type la SNRs has a low temperature (~5000 K: 
Flamilton & Fesen 1988), the Mach number is estimated to be a 
few thousand. Previously, the nature of collisionless heating at 
SNR blast waves has been studied via optical spectroscopy of 
B aimer-dominated shocks (e.g., Raymond et al. 1983; Laming 
et al. 1996; Ghavamian et al. 2001, 2007; Rakowski et al. 
2003; Helder et al. 201 1). However, temperature measurements 
using optical spectra have been subject to large uncertainties 
in atomic cross sections, especially for high-velocity shocks 
( v s > 1000 km s _1 ; e.g., Heng 2010). 

For Balmer-dominated shocks, two main scenarios are sug- 
gested as the origin of the electron heating: (1) lower hybrid 
wave heating in a cosmic-ray precursor (e.g.. Laming 2000; 
Ghavamian et al. 2007) and (2) plasma wave heating due to 
Buneman instabilities formed by reflected non-Maxwellian ions 
(e.g., Cargill & Papadopoulos 1988; Matsukiyo 2010). Both sce- 
narios require a quasi-perpendicular shock, which is unlikely for 
RSs, where the magnetic field is expected to be quasi-parallel to 
the fluid flow due to expansion-induced stretch of the field lines. 
Moreover, there is little evidence for relativistic cosmic rays at 
the RS of Tycho's SNR (Warren et al. 2005), which makes sce- 
nario (a) unlikely in this context. An alternative scenario, the 
cross-shock potential (e.g., Balikhin et al. 1993), which has also 
been suggested as the origin of energetic electrons in GRB af- 
terglows (Gedalin et al. 2008; Sironi & Spitkovsky 2011), may 
apply to our case. In this model, charge separation is created 


at the shock front due to the different gyroradii of ions and 
electrons, forming a potential gap where electrons arriving later 
can be accelerated. This requires no specific orientation of the 
background magnetic field and predicts the self-generation of 
small-scale electromagnetic fields (Gedalin et al. 2008). The 
high Mach number for the RS also supports this analogy with 
GRB shocks. Our observation of efficient collisionless heating 
in the unique environment of a SNR RS suggests that these 
shocks may be fundamentally different from the more widely 
studied Balmer-dominated shocks into ISM material. 

We are thankful to Drs. John D. Raymond and Timothy R. 
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